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Abstract 



After a brief history of neutron stars and supernovae recent developments are dis- 
cussed. Based on modern nucleon-nucleon potentials more reliable equations of state 
for dense nuclear matter have been constructed. Furthermore, phase transitions such 
\^J \ as pion, kaon and hyperon condensation, superfluidity and quark matter can occur 

(^ ' in cores of neutron stars. Specifically, the nuclear to quark matter phase transition 

^^ . and its mixed phases with intriguing structures is treated. Rotating neutron stars 

with and without phase transitions are discussed and compared to observed masses, 
radii and glitches. The observations of possible heavy ~ 2Mq neutron stars in X- 
0^ . ray binaries and QPO's require relatively stiff equation of states and restricts strong 

^^ \ phase transitions to occur at very high nuclear densities only. 
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1 INTRODUCTION AND HISTORY OF NEUTRON STARS 



It is a great pleasure to participate in the 6'th CoUoque Cosmologie at the 
Observatoire de Paris and stand on the same historic spot, where my fellow 
countryman Ole R0mer more than two centuries ago first calculated the speed 
of lightQ. I shall therefore start with a brief history of neutron stars and 
subsequently discuss some very important and recent developments in this 
field at the turn of this millenium. The most important discoveries concerning 
neutron stars are briefly listed in Table I. Most are well known except perhaps 
for the most recent ones, which will be discussed in more detail below. 

Table 1 

Chronological list of important developments related to neutron stars 



Year "Observers" 



Discovery 



1054 


Chinese 


1572 


Tycho Brahe 


1604 


Kepler 


1932 


Chadwick 


1932 


Landau 


1934 


Baade & Zwicky 


1935 


Oppenheimer & Volkoff 


1967 


Bell & Hewish 


1969 




1973 


Hulse & Taylor 


1987 


Neutrino detectors 


1995 


Nijmegen data base 


1996 


RXTE 


1997 


BeppoSAX 



record the Crab Supernova 

observes a Supernova 

observes a Supernova 

discovers the neutron 

suggests the existence of neutron stars 

connects supernovae to gravitational 

collapse of stars to neutron stars 

calculate neutron star structures with simple EoS 

discover first pulsar 

pulsars in Crab and Vela Supernova remnants 

discover first binary pulsar 

collects 19 neutrinos from SN-1987A 

compilation of ^5000 NN cross sections 

leads to "modern" NN potentials and EoS 

discover kHz oscillations (QPO) in X-ray binaries 

Gamma Ray Burst with afterglow at z^l 



^ Actually, he measured and (over-)estimated the light time from the sun, ~ 12mm, 
but did not bother to calculate the speed of light because the very distance to the 
sun was not well known at that time. 



2 HEAVY NEUTRON STARS IN X-RAY BINARIES 



The best determined neutron star masses are found in binary pulsars and all 
lie in the range 1.35 ± O.O4M0 (see Thorsett and Chakrabarty 1999). These 
masses have been accurately determined from variations in their radio pulses 
due to doppler shifts as well periastron advances of their close elliptic orbits 
that are strongly affected by general relativistic effects. One exception is the 
nonrelativistic pulsar PSR J1012+5307 of massQ M = (2.1 ± 0.8)Mq (van 
Paradijs 1998). 

Several X-ray binary masses have been measured of which the heaviest are Vela 
X-1 with M = (1.9 ± 0.2)Mq (Barziv et al, 1999) and Cygnus X-2 with M = 
(1.8 ± OA)Mq (Orosz & Kuulkers 1999). Their Kepler orbits are determined 
by measuring doppler shift of both the X-ray binary and its companion. To 
complete the mass determination one needs the orbital inclination which is 
determined by eclipse durations, optical light curves, or polarization variations 
(see, e.g., van Paradijs). 

The recent discovery of high-frequency brightness oscillations in low-mass X- 
ray binaries provides a promising new method for determining masses and 
radii of neutron stars (see Miller, Lamb, & Psaltis 1998). The kilohertz quasi- 
periodic oscillations (QPO) occur in pairs and are most likely the orbital 
frequencies 



i^QPO = {l/2n) JGM/Rl 



orb ' 



of accreting matter in Keplerian orbits around neutron stars of mass M and 
its beat frequency with the neutron star spin, uqpo — i^s- According to Zhang, 
Strohmayer, & Swank 1997, Kaaret, Ford, & Chen (1997) the accretion can 
for a few QPO's be tracked to its innermost stable orbit, 

R,ns = 6GM/c^ . (2) 

For slowly rotating stars the resulting mass is from Eqs. (1,2) 

kHz 

M ~ 2.2M0 . (3) 

l^QPO 

For example, the maximum frequency of 1060 Hz upper QPO observed in 
4U 1820-30 gives M ~ 2.25M0 after correcting for the Vs — 275 Hz neu- 
tron star rotation frequency. If the maximum QPO frequencies of 4U 1608-52 
{^QPO = 1125 Hz) and 4U 1636-536 {vqpo = 1228 Hz) also correspond to 
innermost stable orbits the corresponding masses are 2.IM0 and 1.9M0. Evi- 



The uncertainties are all 95% conf. limits or ^ 2a 



dence for the innermost stable orbit has only been found for 4U 1820-30, where 
uqpo display a distinct saturation with accretion rate indicating that orbital 
frequency cannot exceed that of the innermost stable orbit. More observations 
like that are needed before a firm conclusion can be made. 

Such large neutron star masses of order ~ 2Mq severely restrict the equation 
of state (EoS) for dense matter as addressed in the following. 



3 MODERN NUCLEAR EQUATION OF STATES 



Recent models for the nucleon-nucleon (NN) interaction, based on the compi- 
lation of more than 5000 NN cross sections in the Nijmegen data bank, have 
reduced the uncertainty in NN potentials. Including many-body effects, three- 
body forces, relativistic effects, etc., the nuclear EoS have been constructed 
with reduced uncertainty allowing for more reliable calculations of neutron 
star properties, see Akmal, Pandharipande, & Ravenhall (1998) and Engvik 
et al. (1997). Likewise, recent realistic effective interactions for nuclear matter 
obeying causality at high densities, constrain the EoS severely and thus also 
the maximum masses of neutron stars, see Akmal, Pandharipande, & Raven- 
hall (1998) and Kalogera & Baym (1996). We have in [14] elaborated on these 
analyses by incorporating causality smoothly in the EoS for nuclear matter 
allowing for first and second order phase transitions to, e.g., quark matter. 

For the discussion of the gross properties of neutron stars we will use the 
optimal EoS of Akmal, Pandharipande, & Ravenhall (1998) (specifically the 
Argonne V18+Sv+ UIX* model- hereafter APR98), which is based on the most 
recent models for the nucleon-nucleon interaction, see Engvik et al. (1997) for 
a discussion of these models, and with the inclusion of a parametrized three- 
body force and relativistic boost corrections. The EoS for nuclear matter is 
thus known to some accuracy for densities up to a few times nuclear saturation 
density no = 0.16 fm~^. We parametrize the APR98 EoS by a simple form for 
the compressional and symmetry energies that gives a good fit around nuclear 
saturation densities and smoothly incorporates causality at high densities such 
that the sound speed approaches the speed of light. This requires that the 
compressional part of the energy per nucleon is quadratic in nuclear density 
with a minimum at saturation but linear at high densities 



£ = Ecomp{n) + S{n){l- 2x) 



2 



u — 2 — s 

= Sou- + Sou^il - 2xf. (4) 

1 + su 

Here, n = Up + Un is the total baryon density, x = np/n the proton fraction 
and u = n/riQ is the ratio of the baryon density to nuclear saturation density. 




The compressional term is in Eq. (4) parametrized by a simple form which 
reproduces the saturation density and the binding energy per nucleon £q = 
15.8MeV at Uq of APR98. The "softness" parameter s ~ 0.2, which gave 
the best fit to the data of APR98 (see Heiselberg & Hjorth- Jensen 1999) is 
determined by fitting the energy per nucleon of APR98 up to densities of 
n ~ AriQ. For the symmetry energy term we obtain 5*0 = 32 MeV and 7 = 0.6 
for the best fit. The proton fraction is given by /3-equilibrium at a given density. 

The one unknown parameter s expresses the uncertainty in the EoS at high 
density and we shall vary this parameter within the allowed limits in the 
following with and without phase transitions to calculate mass, radius and 
density relations for neutron stars. The "softness" parameter s is related to 
the incompressibility of nuclear matter as Kq = 18£^o/(l + s) ~ 200MeV. It 
agrees with the poorly known experimental value (Blaizot, Berger, Decharge, 
& Girod 1995), Kq ~ 180 — 250MeV which does not restrict it very well. From 
{vs/cY = dP/d{nS), where P is the pressure, and the EoS of Eq. (4), the 
causality condition Cg < c requires 



0.13, (5) 



where rrin is the mass of the nucleon. With this condition we have a causal 
EoS that reproduces the data of APR98 at densities up to 0.6 ~ 0.7 fm~^. In 
contrast, the EoS of APR98 becomes superluminal at n ~ 1.1 fm~^. For larger 
s values the EoS is softer which eventually leads to smaller maximum masses 
of neutron stars. The observed M ~ 1.4M0 in binary pulsars restricts s to be 
less than 0.4 — 0.5 depending on rotation as shown in calculations of neutron 
stars below. 

In Fig. 1 we plot the sound speed {vs/cf for various values of s and that 
resulting from the microscopic calculation of APR98 for /9-stable pn-matter. 
The form of Eq. (4), with the inclusion of the parameter s, provides a smooth 
extrapolation from small to large densities such that the sound speed Vg ap- 
proaches the speed of light. For s = 0.0 (s = 0.1) the EoS becomes superlu- 
minal at densities of the order of 1 (6) fm^^. 

The sound speed of Kalogera & Baym (1996) is also plotted in Fig. 1. It jumps 
discontinuously to the speed of light at a chosen density. With this prescription 
they were able to obtain an optimum upper bound for neutron star masses and 
obey causality. This prescription was also employed by APR98, see Rhoades 
& Rufiini (1974) for further details. The EoS is thus discontinuously stiffened 
by taking f ^ = c at densities above a certain value ric which, however, is lower 
than Hs = 5no where their nuclear EoS becomes superluminal. This approach 
stiffens the nuclear EoS for densities Uc < n < Ug but softens it at higher den- 
sities. Their resulting maximum masses lie in the range 2.2M0^M^2.9M0. 
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Fig. 1. {vs/cf for /3-stable pn-matter for s = 0,0.1,0.2,0.3, the results of APR98, 
and for the patched EoS of Kalogera & Baym (1996) which shows a discontinuous 

(Vs/C)'. 

Our approach however, incorporates causality by reducing the sound speed 
smoothly towards the speed of light at high densities. Therefore our approach 
will not yield an absolute upper bound on the maximum mass of a neutron star 
but gives reasonable estimates based on modern EoS around nuclear matter 
densities, causality constraints at high densities and a smooth extrapolation 
between these two limits (see Fig. 1). 



At very high densities particles are expected to be relativistic and the sound 
speed should be smaller than the speed of light, v'j ^ c^/3. Consequently, 
the EoS should be even softer at high densities and the maximum masses we 
obtain with the EoS of (4) are likely to be too high estimates. 



4 PHASE TRANSITIONS 



The physical state of matter in the interiors of neutron stars at densities above 
a few times normal nuclear matter densities is essentially unknown and many 
first and second order phase transitions have been speculated upon. We will 
specifically study the hadron to quark matter transition at high densities, 
but note that other transitions as, e.g., kaon and/or pion condensation or the 
presence of other baryons like hyperons also soften the EoS and thus further 
aggravate the resulting reduction in maximum masses. Hyperons appear at 
densities typically of the order 2no and result in a considerable softening of the 
EoS, see e.g., Balberg, Lichenstadt, & Cook (1998). Typically, most equations 
of state with hyperons yield masses around 1.4 — 1.6Mq. Here however, in 
order to focus on the role played by phase transitions in neutron star matter, 
we will assume that a phase transition from nucleonic to quark matter takes 
place at a certain density. Since we do not have a fully reliable theory for 
the quark matter phase, we will for simplicity employ the bag model in our 
actual studies of quark phases and neutron star properties. In the bag model 
the quarks are assumed to be confined to a finite region of space, the so-called 
'bag', by a vacuum pressure B. Adding the Fermi pressure and interactions 
computed to order ag = g^/Air, where g is the QCD coupling constant, the 
total pressure for 3 massless quarks of flavor f = u,d, s, is (see Kapusta (1988), 



47r^ TT 



a^) - 5 + Pe + Pm . (6) 



where Pe,/i are the electron and muon pressure, e.g., Pe = /^e/^^TT^. A Fermi 
gas of quarks of flavor i has density Uj = k^^ji^^ ^ due to the three color states. 
A flnite strange quark mass have minor effect on the EoS since quark chemical 
potentials ^q^mj^/?) typically are much larger. The value of the bag constant 
B is poorly known, and we present results using two representative values, 
B = 150 MeVfm-3 and B = 200 MeVfm'^. We take a^ = 0.4. However, 
similar results can be obtained with smaller ag and larger B (Madsen 1998). 

The quark and nuclear matter mixed phase described in Glendenning (1992) 
has continuous pressures and densities due to the general Gibbs criteria for 
two-component systems. There are no flrst order phase transitions but at most 
two second order phase transitions. Namely, at a lower density, where quark 
matter flrst appears in nuclear matter, and at a very high density (if gravita- 
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Fig. 2. The quark and nuclear matter structure in a quarter of a typical IAMq solar 
mass neutron star. The typical sizes of structures are a few Fermi's but have been 
scaled up by about 16 orders of magnitudes to be seen. 

tionally stable), where all nucleons are finally dissolved into quark matter. This 
mixed phase does, however, not include local surface and Coulomb energies 
of the quark and nuclear matter structures. If the interface tension between 
quark and nuclear matter is too large, the mixed phase is not favored ener- 
getically due to surface and Coulomb energies associated with forming these 
structures (Heiselberg, Pethick, & Staubo 1993). The neutron star will then 
have a core of pure quark matter with a mantle of nuclear matter surround- 
ing it and the two phases are coexisting by a first order phase transition or 
Maxwell construction. For a small or moderate interface tension the quarks 
are confined in droplet, rod- and plate-like structures as found in the inner 
crust of neutron stars (Lorenz, Ravenhall, & Pethick 1993). 



5 MASSES AND RADII OF NEUTRON STARS 



In order to obtain the mass and radius of a neutron star, we have solved 
the Tolman-Oppenheimer-Volkov equation with and without rotational correc- 
tions following the approach of Hartle (1967). The equations of state employed 
are given by the pn-matter EoS with s = 0.13, 0.2, 0.3, 0.4 with nucleonic de- 



grees of freedom only. In addition we have selected two representative values for 
the bag-model parameter B, namely 150 and 200 MeVfm"^ for our discussion 
on eventual phase transitions. The quark phase is linked with our pn-matter 
EoS from Eq. (4) with s = 0.2 through either a mixed phase construction or 
a Maxwell construction, see Heiselberg and Hjorth- Jensen (1999) for further 
details. For B = 150 MeVfm"^, the mixed phase begins at 0.51 fm~^ and the 
pure quark matter phase begins at 1.89 fm~^. Finally, for B = 200 MeVfm"^, 
the mixed phase starts at 0.72 fm~^ while the pure quark phase starts at 2.11 
fm~^. In case of a Maxwell construction, in order to link the pn and the quark 
matter EoS, we obtain for B = 150 MeVfm"^ that the pure pn phase ends 
at 0.92 fm~^ and that the pure quark phase starts at 1.215 fm~^, while the 
corresponding numbers for B = 200 MeVfm"^ are 1.04 and 1.57 fm^^. 

None of the equations of state from either the pure pn phase or with a mixed 
phase or Maxwell construction with quark degrees of freedom, result in stable 
configurations for densities above ~ lOno, implying thereby that none of the 
stars have cores with a pure quark phase. The EoS with pn degrees of freedom 
have masses M^2.2M0 when rotational corrections are accounted for. With 
the inclusion of the mixed phase, the total mass is reduced since the EoS 
is softer. However, there is the possibility of making very heavy quark stars 
for very small bag constants. For pure quark stars there is only one energy 
scale namely B which provides a homology transformation (Madsen 1998) 
and the maximum mass is Mmax = 2.OM0(58MeVfm~^/_B)^/^ (for as = 0). 
However, for i?^58MeVfm~^ a nuclear matter mantle has to be added and for 
-B^58MeVfm~^ quark matter has lower energy per baryon than ^^Fe and is 
thus the ground state of strongly interacting matter. Unless the latter is the 
case, we can thus exclude the existence of 2.2 — 2.3Mq quark stars. 

In Fig. 5 we show the mass-radius relations for the various equations of state. 
The shaded area represents the allowed masses and radii for uqpq = 1060 Hz 
of 4U 1820-30. Generally, 



2GM<R<{j-^-^ , (7) 

^QPO/ 



47rV; 



where the lower limit insures that the star is not a black hole, and the upper 
limit that the accreting matter orbits outside the star, R < Rort- Furthermore, 
for the matter to be outside the innermost stable orbit, R > Rms = QGM, 
requires that 

M < ^ + °-^'^" ^ 2.2M,(1 + 0.75j)>''' 



where j = 2'kcVsI /M"^ is a dimensionless measure of the angular momentum 
of the star with moment of inertia /. The upper limit in Eq. (8) is the mass 
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Fig. 3. Total mass M as function of central density n^ for various values of s (upper 
panel) and the bag parameter B (lower panel) for both a mixed phase and a Maxwell 
constructed EoS with s = 0.2 in Eq. 4). In addition we include also the rotational 
corrections for the pure pn-case with s = 0.2 and the mixed phase construction 
for B = 200 MeVfm"^. For the Maxwell construction which exhibits a first order 
phase transition, in the density regions where the two phases coexist, the pressure is 
constant, a fact reflected in the constant value of the neutron star mass. All results 
are for /^-stable matter. Note also that for the upper panel, the EoS for s = 0.3 and 
s = 0.4 start to differ from those with s = 0.13, 0.2 at densities below 0.2 fm~^. 
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when vqpo corresponds to the innermost stable orbit. According to Zhang, 
Smale, Strohmayer & Swank (1998) this is the case for 4U 1820-30 since vqpo 
saturates at ~ 1060 Hz with increasing count rate. The corresponding neutron 
star mass is M ~ 2.2 — 2.3M0 which leads to several interesting conclusions 
as seen in Fig. 5. Firstly, the stiffest EoS allowed by causality (s ~ 0.13 — 0.2) 
is needed. Secondly, rotation must be included which increase the maximum 
mass and corresponding radii by 10-15% for Vs ~ 300 Hz. Thirdly, a phase 
transition to quark matter below densities of order ~ Sng can be excluded, 
corresponding to a restriction on the bag constant 5^200 MeVfm"^. 

These maximum masses are smaller than those of APR98 and Kalogera & 
Baym (1996) who, as discussed above, obtain upper bounds on the mass of 
neutron stars by discontinuously setting the sound speed to equal the speed 
of light above a certain density, Hf.. By varying the density rZc = 2 — > 5no the 
maximum mass drops from 2.9 -^ 2.2Mq. In our case, incorporating causality 
smoothly by introducing the parameter s in Eq. (4), the EoS is softened at 
higher densities in order to obey causality, and yields a maximum mass which 
instead is slightly lower than the 2.2M0 derived in APR98 for nonrotating 
stars. 

If the QPOs are not from the innermost stable orbits and one finds that 
even accreting neutron stars have small masses, say like the binary pulsars, 
M^1.4M0, this may indicate that heavier neutron stars are not stable. There- 
fore, the EoS is soft at high densities s^O.4 or that a phase transition occurs 
at a few times nuclear matter densities. For the nuclear to quark matter tran- 
sition this would require B < 80 MeVfm"^ for s = 0.2. For such small bag 
parameters there is an appreciable quark and nuclear matter mixed phase in 
the neutron star interior but even in these extreme cases a pure quark matter 
core is not obtained for stable neutron star configurations. 



6 GLITCHES 



Younger pulsars rotate and slow down rapidly. Some display sudden speed 
ups referred to as glitches. We shall first discuss such standard glitches and 
subsequently continue to giant glitches and other characteristic glitch behavior 
for spin frequencies close to the critical ones, where first order phase transitions 
occur at densities present right at the core of neutron stars. 
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Fig. 4. Neutron star masses vs. radius for the EoS of Eq. (1) with softness 
s=0. 13,0.2,0.3,04, with increasing values of s from top to bottom for the full curves. 
Phase transitions decrease the maximum mass whereas rotation increases it. The 
shaded area represents the neutron star radii and masses allowed (see text and Eqs. 
1-3) for orbital QPO frequencies 1060 Hz of 4U 1820-30. 

6. 1 Core quakes and glitches 



The glitches observed in the Crab, Vela, and a few other pulsars are probably 
due to quakes occurring in solid structures such as the crust, superfluid vortices 
or possibly the quark matter lattice in the core [16]. As the rotating neutron 
star gradually slows down and becomes less deformed, the rigid component 
is strained and eventually cracks/quakes and changes its structure towards 
being more spherical. 



The moment of inertia of the rigid component, Ic-, decreases abruptly and its 
rotation and pulsar frequency increases due to angular momentum conserva- 
tion resulting in a glitch. The observed glitches are very small AVL/VL ~ 10~*. 
The two components slowly relaxate to a common rotational frequency on a 
time scale of days (healing time) due to superfluidity of the other component 
(the neutron liquid). The healing parameter Q = Id hot measured in glitches 
reveals that for the Vela and Crab pulsar about ~3% and ~96% of the moment 
of inertia is in the rigid component respectively. 
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If the crust were the only rigid component the Vela neutron star should be 
almost all crust. This would require that the Vela is a very light neutron star - 
much smaller than the observed ones which all are compatible with ~ IAMq. 
If we by the lattice component include not only the solid crust but also the 
protons in nuclear matter (NM) (which is locked to the crust due to magnetic 
fields), superfluid vortices pinned to the crust [27] and the solid QM mixed 
phase 

Ic = hrust + Ip + Isv + IqM i (9) 

we can better explain the large Ic for the Crab. The moment of inertia of 
the mixed phase is sensitive to the EoS's used. For example, for a quadratic 
NM EoS [16] decreasing the Bag constant from 110 to 95 MeVfm^^ increases 
Iclhotai from ~ 20% to ~ 70% for a IAMq neutron star - not including possible 
vortex pinning. The structures in the mixed phase would exhibit anisotropic 
elastic properties, being rigid to some shear strains but not others in much 
the same way as liquid crystals. Therefore the whole mixed phase might not 
be rigid. 

The energy released in glitches every few years are too large to be stored in the 
crust only. The recurrence time for large quakes, tc, is inversely proportional 
to the strain energy [27], which again is proportional to the lattice density and 
the Coulomb energy 

t.-^oc- . 10 

a^ a 

Since the lattice distance a is smaller for the quark matter droplets and their 
charge larger than for atoms in the crust, the recurrence time is shorter in 
better agreement with measurements of large glitches. 

Detecting core and crust quakes separately or other signs of three components 
in glitches, indicating the existence of a crust, superfluid neutrons and a solid 
core, would support the idea of the mixed quark and nuclear matter mixed 
phase. However, magnetic fleld attenuation is expected to be small in neutron 
stars and therefore magnetic flelds penetrate through the core. Thus the crust 
and core lattices as well as the proton liquid should be strongly coupled and 
glitch simultaneously. 



6.2 Backbending and giant glitches 



In [11] the moment of inertia is found to "backbend" as function of angu- 
lar velocity. The moment of inertia of some deformed nuclei [24,17] may also 
backbend when the coriolis force exceeds the pairing force breaking the pair- 
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ing whereby the nucleus reverts from partial superfluidity to a rigid rotor. 
However, in the limit of large nuclear mass number such backbending would 
disappear. Instead pairing may lead to superfluidity in bulk [27]. A backbend- 
ing phenomenon in neutron stars, that appears to be similar to backbending 
in nuclei, can occur in neutron stars although the physics behind is entirely 
different. If we soften the EoS significantly at a density near the central den- 
sity of the neutron star, a non-rotating neutron star can have most of its core 
at high densities where the soft EoS determines the profile. A rapidly rotating 
star may instead have lower central densities only probing the hard part of the 
EoS. Thus the star may at a certain angular velocity revert from the dense 
phase to a more dilute one and at the same time change its structure and 
moment of inertia discontinuosly. Such a drastic change in moment of inertia 
at some angular velocity will cause a giant glitch as found in [11]. 

The phenomenon of neutron star backbending is related to the double maxi- 
mum mass for a neutron and quark star respectively as discussed in Ref. [14]. 
Both phenomena occur when the "softening" or density discontinuity exceeds 

3 

£2~-^i =^ Two maximum masses and Giant glitch when u; ~ cjq- (H) 

Similarly, there is also a discontinuity in the moment of inertia when p2/Pi~3/2. 
The neutron star may continue to slow down in its unstable structure, i.e. 
"super-rotate" , before reverting to its stable configuration with a dense core. 
As for the instabilities of Eq. (11) this condition changes slightly for a more 
general EoS and when general relativity is included. Neutron stars with a 
mixed phase do not have a first order phase transition but may soften their 
EoS sufficiently that a similar phenomenon occurs [11]. We emphasize, how- 
ever, that the discontinuous jump in moment of inertia is due to the drastic 
and sudden softening of the EoS near the central density of neutron stars. It 
is not important whether it is a phase transition or another phenomenon that 
causes the softening. 



6.3 Phase transitions in rotating neutron stars 



During the last years, and as discussed in the previous sections, interesting 
phase transitions in nuclear matter to quark matter, mixed phases of quark 
and nuclear matter [10,16], kaon [20] or pion condensates [1], neutron and pro- 
ton superfluidity [8], hyperonic matter [33,32,21] crystalline nuclear matter, 
magnetized matter, etc., have been considered. Here we consider the inter- 
esting phenomenon of how the star and in particular its moment of inertia 
behaves near the critical angular velocity, Qo, where the core pressure just 
exceeds that needed to make a phase transition. 
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As described in detail in [14] the general relativistic equations for slowly rotat- 
ing stars can be solved even with first order phase transitions. The resulting 
moment of inertia have the characteristic behavior for i^^fio (see also Fig. 
(5)) 

-|c.(l-|p + ...). (12) 

For the two incompressible fluids with momentum of inertia, the small expan- 
sion parameters are ci = cUq and C2 = {h/2)ujl{e2/si — l)/(3 — 262/^1 — ^0)'^^^) 
for i7 > r^o the C2 term is absent. For a Bethe- Johnson polytropic EoS we find 
from Fig. (5) that C2 — 0.07 ~ 2.2uJq. Generally, we find that the coefficient 
C2 is proportional to the density difference between the two coexisting phases 
and to the critical angular velocity to the third power, C2 ~ (^2/^1 ~ l)^o- The 
scaled critical angular velocity ooq can at most reach unity for submillisecond 
pulsars. 

To make contact with observation we consider the temporal behavior of an- 
gular velocities of pulsars. The pulsars slow down at a rate given by the loss 
of rotational energy which we shall assume is proportional to the rotational 
angular velocity to some power (for dipole radiation n = 3) 

With the moment of inertia given by Eq. (12) the angular velocity will then 
decrease with time as 



ClT^ - C2 



*' jq \ a Zq y ^ 'n 



1 /. L ,h 



for t > to. Here, the time after formation of the pulsar is, using Eq. (13), related 
to the angular velocity as t ~ to(i^o/^)"~"'^ and to = lo/Hn — 1)C0.q~^) for 
n > 1, is the critical time where a phase transition occurs in the center. For 
earlier times t < to there is no dense core and Eq. (14) applies when setting 
C2 = The critical angular velocity is Qo = ujoy/^irei ~ 6kHz for the Bethe- 
Johnson EoS [4], i.e. comparable to a millisecond binary pulsar. Applying 
these numbers to, for example, the Crab pulsar we find that it would have 
been spinning with critical angular velocity approximately a decade after the 
Crab supernova explosion, i.e. to ~ 10 years for the Crab. Generally, to oc f2o~" 
and the timescale for the transients in Cl as given by Eq. (13) may be months 
or centuries. In any case it would not require continuous monitoring which 
would help a dedicated observational program. 
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The braking index depends on the second derivative /" = dI/(PQ of the 
moment of inertia and thus diverges (see Fig. (1)) as Q approaches Qq from 
below 

For r2 > r^o the term with C2 is absent. The observational braking index n{Q) 
should be distinguished from the theoretical exponent n appearing in Eq. (12). 
Although the results in Eqs. (13) and (14) were derived for the pulsar slow 
down assumed in Eq. (12) both Cl and n{Q) will generally display the y/t — to 
behavior for t > to as long as the rotational energy loss is a smooth function 
of fi. The singular behavior will, however, be smeared on the pulsar glitch 
"healing" time which in the case of the Crab pulsar is of order weeks only. 

We now discuss possible phase transitions in interiors of neutron stars. The 
quark and nuclear matter mixed phase described in [10] has continuous pres- 
sures and densities. There are no first order phase transitions but at most 
two second order phase transitions. Namely, at a lower density, where quark 
matter first appears in nuclear matter, and at a very high density (if gravita- 
tionally stable), where all nucleons are finally dissolved into quark matter. In 
second-order phase transitions the pressure is a continuous function of density 
and we find a continuous braking index. This mixed phase does, however, not 
include local surface and Coulomb energies of the quark and nuclear matter 
structures. As shown in [16] there can be an appreciable surface and Coulomb 
energy associated with forming these structures and if the interface tension 
between quark and nuclear matter is too large, the mixed phase is not favored 
energetically. The neutron star will then have a core of pure quark matter with 
a mantle of nuclear matter surrounding it and the two phases are coexisting 
by a first order phase transition. For a small or moderate interface tension 
the quarks are confined in droplet, rod- and plate-like structures as found in 
the inner crust of neutron stars [22]. Due to the finite Coulomb and surface 
energies associated with forming these structures, the transitions change from 
second order to first order at each topological change in structure. If a Kaon 
condensate appears it may also have such structures [12]. Pion condensates 
[7], crystalline nuclear matter [1], hyperonic or magnetized matter, etc. may 
provide other first order phase transitions. 

Thus, if a first order phase transitions is present at central densities of neutron 
stars, it will show up in moments of inertia and consequently also in angular 
velocities in a characteristic way. For example, the slow down of the angular 



velocity has a characteristic behavior f2 ~ C2a/1 — t/to and the braking index 

diverges as n{il) ~ C2/J1 — H'^/Hq (see Eqs. (14) and (15)). The magnitude of 
the signal generally depends on the density difference between the two phases 
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and the critical angular velocity c^o = ^o/^/2rcei such that C2 ~ (£2/^1 — l)^o- 
The observational consequences depend very much on the critical angular 
velocity Dq, which depends on the EoS employed, at which density the phase 
transition occurs and the mass of the neutron star. 

We encourage a dedicated search for the characteristic transients discussed 
above. As the pulsar slows down over a million years, its central densities 
spans a wide range of order Itiq (see Fig. 5). As we are interested in time 
scales of years, we must instead study all ~ 1000 pulsars available. By study- 
ing the corresponding range of angular velocities for the sample of different 
star masses, the chance for encountering a critical angular velocity increases. 
Eventually, one may be able to cover the full range of central densities and find 
all first order phase transitions up to a certain size determined by the experi- 
mental resolution. Since the size of the signal scales with VL^ the transition may 
be best observed in rapidly rotating pulsars such as binary pulsars or pulsars 
recently formed in supernova explosion and which are rapidly slowing down. 
Carefully monitoring such pulsars may reveal the characteristic behavior of 
the angular velocity or braking index as described above which is a signal of 
a first order phase transition in dense matter. 



7 SUMMARY 



Modern nucleon-nucleon potentials have reduced the uncertainties in the cal- 
culated EoS. Using the most recent realistic effective interactions for nuclear 
matter of APR98 with a smooth extrapolation to high densities including 
causality, the EoS could be constrained by a "softness" parameter s which 
parametrizes the unknown stiffness of the EoS at high densities. Maximum 
masses have subsequently been calculated for rotating neutron stars with and 
without first and second order phase transitions to, e.g., quark matter at high 
densities. 

The calculated bounds for maximum masses leaves two natural options when 
compared to the observed neutron star masses: 

• Case I: The large masses of the neutron stars in QPO 4U 1820-30 (M = 
2.3Mq), PSR J1012+5307 (M = 2.1 ± OAMq), Vela X-1 (M = 1.9 ± 
O.IMq), and Cygnus X-2 (M = 1.8 ± 0.2Mq), are confirmed and comple- 
mented by other neutron stars with masses around ~ 2Mq. 

As a consequence, the EoS of dense nuclear matter is severely restricted 
and only the stiffest EoS consistent with causality are allowed, i.e., softness 
parameter 0.13 < s^O.2. Furthermore, any significant phase transition at 
densities below < 5no can be excluded. 

That the radio binary pulsars all have masses around IAMq is then proba- 
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Fig. 5. Central density (in units of Eq), radii of the neutron star R and its dense 
core Rq, moment of inertia, its derivative I' /I = dl/du'^ /I and the braking index 
are shown as function of the scaled angular velocity w^ = f2^/(27rei). The rotating 
neutron star has mass IAMq and a Bethe- Johnson like polytropic EoS with a first 
order phase transition taking place at density ei = 3.2eo to £2 = 4eo- 

bly due to the formation mechanism in supernovae where the Chandrasekhar 
mass for iron cores are ~ 1.5Mq. Neutron stars in binaries can subsequently 
acquire larger masses by accretion as X-ray binaries. 
• Case II: The heavy neutron stars prove erroneous by more detailed obser- 
vations and only masses like those of binary pulsars are found. 

If accretion does not produce neutron stars heavier than ^IAMq, this 
indicates that heavier neutron stars simply are not stable which in turn 
implies a soft EoS, either s > 0.4 or a significant phase transition must 
occur already at a few times nuclear saturation densities. 



Surface temperatures can be estimated from spectra and from the measured 
fluxes and known distances, one can extract the surface area of the emitting 
spot. This gives unfortunately only a lower limit on the neutron star size, R. 
If it becomes possible to measure both mass and radii of neutron stars, one 



can plot an observational (M, R) curve in Fig. (5), which uniquely determines 
the EoS for strongly interacting matter at zero temperature. 

Pulsar rotation frequencies and glitches are other promising signals that could 
reveal phase transitions. Besides the standard glitches also giant glitches were 
mentioned and in particular the characteristic behavior of angular velocities 
when a first order phase transition occurs right in the center of the star. 

It is, unfortunately, impossible to cover all the interesting and recent devel- 
opments concerning neutron stars in these proceedings so for more details we 
refer to [14] and Refs. therein. 

Acknowledgements: to my collaborators M. Hjorth- Jensen and C.J. Pethick 
as well as G. Baym, F. Lamb and V. Pandharipande. 
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